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This paper examines the stability of a previously proposed version of the shear-transformation- 
zone (STZ) theory of plasticity where the total STZ population is determined by an effective tem- 
perature and compares it to experimental results for a metallic glass. In particular, the addition 
of effective temperature dynamics to the shear transformation zone theory leads to the existence of 
a range of strain rates for which the strain localizes into shear bands. Yet while the steady-state 
results qualitatively agree, the instability of the dynamics of the system while loading begins better 
describes the experimental observations. 
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I. INTRODUCTION 

Recently, Langer[l] derived a finite-temperature ver- 
sion of shear-transformation-zonc (STZ) theory wherein 
the equilibrium STZ density was determined by an ef- 
fective temperature and qualitatively agreed with that 
experimental results on metallic glasses done by Lu, et 
al.[2] The rheology of both is characterized by viscos- 
ity that is Newtonian at low strain rates and decreases 
as the strain rate increases. In the theory, the Newto- 
nian viscosity at low strain rates can be explained by 
thermal creep, while the decrease in viscosity at higher 
strain rates results from the stress-driven dynamics of the 
zones. The experiments end at high strain rates at which 
the material exhibits shear localization during deforma- 
tion. Towards the end of the paper, Langer hypothesized 
that under certain conditions, the equations could repro- 
duce this shear banding. This paper is an analysis of 
those shear-localization properties. 

Stressed metallic glasses typically fail by shear local- 
ization but the underlying physical cause is under debate. 
While the theories of adiabatic shear bands, in which the 
shear banding instability is caused by the deformation- 
induced increase in the thermal temperature, have been 
successfully applied to crystalline metals [3], experimen- 
tal reports on metallic glasses are in conflict as to 
whether the temperature increases produced by defor- 
mation are sufficiently high enough to reach the melting 
point. [4, 5, 6] 

It has been hypothesized that the configurational de- 
grees of freedom in a non-equilibrium system, such as a 
metallic glass, can be described by an effective temper- 
ature that does not equal the temperature of the heat 
bath. Plastic rearrangements, such as those during de- 
formation, could increase this effective temperature up 
to a value near an effective melt temperature, and that 
could be a mechanism that induces shear localization. 
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despite the nominal temperature of the material staying 
below the melting point. This will be shown to be the 
case for the STZ model described in [1]. 

Such an effective temperature would play a similar role 
as the free volume in other shear transformation zone 
and flow-defect models, such as Spaepen[7] — a state vari- 
able describing the disorder of a system — but need not 
be physically identical. Recent simulations by Falk and 
Shi [8, 9] have suggested that fluctuations of the free vol- 
ume may be be more important than its absolute value 
in differentiating regions of shear localization from the 
bulk. The effective temperature in the STZ model is a 
way of parametrizing these fluctuations. 

In addition to steady-state shear bands, as discussed 
in [1], Huang et al. showed a shear band can form during 
the initial transients of a controUed-strain-rate experi- 
ment, even when the only possible steady-state solution 
is uniform. While they examined Spaepen free-volume 
model, a similar examination of the STZ model with ef- 
fective temperature will also show instabilities due to the 
initial response. This type of instability will be shown to 
better explain the experimental results than a steady- 
state one. 

This paper begins by quantifying the conditions for 
steady-state shear localization in the shear transforma- 
tion zone model as defined by Langer in [1] and then 
repeats the comparison to the experiments done on the 
metallic glass Zr4i.2Tii3.8Cui2.5NiioBe22.5 (Vitreloy 1) 
by Lu, et al.[2] with particular attention to the localiza- 
tion properties of the material. First, I review the shear 
transformation zone model as described in [1]. In Sec- 
tion III, I examine the steady-state linear stability of the 
model and make some simple comparisons to the some of 
the metallic glass data. The nonlinear contributions to 
the stability and the dynamics of the initial loading are 
included in Section IV. Finally, in Section V, the results 
are summarized and suggestions for future research are 
made. 
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II. STZ EQUATIONS 



where Eq = AePrioc and 



STZ theory hypothesizes that the plastic strain of a 
material under stress is localized in small regions. In 
these regions, called shear transformation zones, the lo- 
cal configuration of the material can be modeled as a 
two state system. Under no stress, each configuration is 
equally stable, but when stress is applied, one orienta- 
tion is preferred over the other. The transition of a STZ 
from one state to the other, called a "flip", results in a 
plastic strain increase Ae of order 1. This increment of 
plastic strain is assumed to be purely deviatoric and non- 
dilational. Once a zone flips, it is incapable of flipping 
any further in that direction, an important distinction 
from previous models. 

Consider a two-dimensional system to which is applied 
a deviatoric stress a. The deviatoric stress tensor can 
be described by a magnitude given by the second stress 

invariant a = ■^aijO'ij . Denoting the density of zones 

that are oriented with (against) a positive stress as 
(ri_), the plastic strain rate can be written as 



Aef 



{R{+s)n- - i?(-s)n+) . 



(1) 



where tq is the characteristic time scale for flipping, on 
the order of the atomic vibrational frequency R{±s) de- 
scribes the stress dependence of the rates of flipping 
=F — *■ ±, and P is the area of an STZ, on the order of a 
couple of atoms. The stress s = ^, where p, is, in general, 
a characteristic of the rates and will be, for the particular 
rates used in this paper, approximately the yield stress. 
The populations of n± are governed by 

Toh± = R{±s)n^ - R{Ts)n± + {T + p) - n±) . (2) 

The first two terms describe the flipping between states 
while the last term describes the creation and annihila- 
tion of zones. The creation and annihilation are driven by 
a rate determined both from the energy dissipated from 
plastic work, r,[10] and from thermal effects, p.[l, If] 
The ability of the total number of STZs to change in 
time is what allows plastic flow at large stresses. 

In the absence of stress, the STZ populations reach 
a steady state value of Acr- In [1], Langer posits that 
this population is described by a Boltzmann distribution 
Acr = ^^e~^^^'^'^'^'ff , governed by an effective temper- 
ature Teff, and a characteristic energy Ez- Ez is the 
energy required for the formation of an STZ, which is ap- 
proximately the cost of a zone due to the the elastic per- 
turbation of the surrounding material, Ez ~ ksTz ~ pl^, 
where p is the shear modulus. 

It is convenient to rewrite the STZ dynamics in terms 
of A = — (nj_ + n_), the total density of STZs, and 

the relative bias of STZs. In those terms, 



To 



-AC(s) (r(s) - to) , 



(3) 



T{s) 



R{+s) - R{-s) 
R{+s) + R{~sy 



The total STZ population A evolves according to 



A 



-{V + p) (e-i/x_A), 

To V / 



(4) 
(5) 

(6) 



where x = "r^' ^^^^ shown in Section II B, when 
the effective temperature dynamics are detailed, that the 
timescale associated with the equilibration of the STZ 
density to e~^^^ is much faster than that associated with 
the effective temperature changes. Ergo for the present 
purposes we shall assume A = e~^/-^ at all times. In this 
limit, the rh equation simplifies to 

TO = — [2C(s)(r(s)-TO)-TO(r + p)]. (7) 

To 

Following the logic of [ ..:], and done in detail in [1], 
the dissipation F can be determined from energy balance 
arguments. Balancing the rate of plastic work with the 
rate of change of the recoverable energy stored in the 
STZs and the dissipation yields 



lit 



2a£P' = ^+^££AF. 



To 



(8) 



Using the fact that the dissipated energy should always 
be positive shows F to have the form 



F(s, to) = 

2C{s) (T{s) - m) [s~T-^{m)) +mpT-^{m) 
1 — mT~^{m) 



(9) 



Since R{—s) is a monotonically increasing function of s, 
T(s) will be also, therefore T^^{m) is both well-defined 
and monotonically increasing. This creates a dynamical 
upper bound on to where ■m'T~^{m) = 1 and the dissi- 
pation becomes infinitely large. Combining Eqs.(7) and 
(9) yields 

1 2C(s) (T(s) — to) (1 — TOs) — mp 
m = : ^ , , — ^ . (10) 



To 



1 — toT ^(to) 



For a system under controlled stress, Eq. (10) is solv- 
able, yielding a steady-state value of to 



rrissis) = — 



1 + sT{s) 



2C{s) 



1 + sT{s) + 



2C{s) 



(11) 
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In the limit p — > 0, where thermal effects are absent, 
Eq. (11) simplifies to the same results as the athermal 
STZ model: a family of jammed solutions where m = 
T(s) and e^' = for \s\ < Sy and a family of STZ flowing 
solutions where m- = 7 and e^' 7^ for |s| > Sy where Sy 
is defined by SyT(sy) = 1. The thermal effects smooth 
this transition by causing slow plastic creep for s < Sy. 



Also the average rate factor C(s) is itself a sharply 
increasing function of stress. It is low for small stresses 
then rapidly increases to ^ when the stress is roughly 
ln(a) . As the temperature decreases, not only does C(0) 
decrease, the transition to C ^ becomes sharper and 
occurs at a higher stress. Eventually, in the limit T ^ 0, 
the total rates will be zero for all finite stress, and such 
a material would exhibit no plasticity. 



Rate factors 



To progress further, wc must get more precise about 
the forms of the rate factors p{T) and R{s). For the 
thermal frequency, p(T) should have the form 



P{T) 



-a(T) 



(12) 



Since at the temperatures of interest, the thermal noise 
frequencies are considerably lower than the atomic vi- 
brational frequency Tq^, a{T) will be greater than 1. 
Determining the actual form of a{T) from a theory is 
beyond the scope of this paper, but to compare to exper- 
iments I shall follow the lead of [] ] and use a Cohen- Grest 
formula[13] as a phcnomenological fit for a{T), 



a{T) 



Tr 



(13) 



For the rate factors i?(s), at stress, any flipping 
should be the result of thermal kicking, and assuming 
the energy barrier for an STZ flip is the same as that in 
creation and annihilation, -R(O) = p{T). Then the sim- 
plest way to add a well-behaved stress dependence is to 
claim that the activation energies in the rates depend ex- 
ponentially on s, ensuring they are always positive. This 
gives R{s) the general form 



R{s) = exp(— a(T) exp(— 7s)) , 



(14) 



For simplicity, I will make the same assumption implied 
in [1] and [12] and use 7=1. Consequently, C and T 
become 

C{s) = exp(— Q!COsh(s))cosh(Q;sinh(s)) (15) 
T(s) = tanh(asinh(s)). (16) 

Remember that the STZ yield stress Sy is defined to sat- 
isfy the condition that SyT{sy) = 1. For a{T) > 2 (which 
for the metallic glasses will turn out to be T < 1000 K), 
Sy ^ jl making /i the yield stress. 

This expression for the transition rate will be adequate 
for the present purposes, but some small issues may arise. 
One possible problem is that these rate factors lead max- 
imum strain rate; as s 00, C{s) i and T ^ 1, limit- 
ing the strain rate to at most ^ exp^— , where Xoo 

is the maximum value of the effective temperature and 
will be described in the next section. 



B. Effective Temperature Dynamics 

As stated earlier, the total STZ population is described 
by an effective temperature %. Because x describes the 
configuration of the system, it increases from plastic 
work, and decreases as the system over time anneals. 

This effective temperature is not simply the free vol- 
ume. In particular, some molecular dynamics simulations 
hint that the free volume may decrease in areas under- 
going plastic strain[9]. Such local free volume decreases 
can still correspond to an increase of the effective tem- 
perature. 

The effective temperature dynamics are determined by 



X = 



To Co 



Ar(s,TO)(xoo-x) + Kpe ^'^ 



T 



g , we will relax that assumption, simply 
is on the order of Tg. 



(17) 

where co is a dimensionless measure of the heat capac- 
ity in a volume the size of an STZ. The first term, the 
heating, is proportional to the total dissipation rate from 
plastic deformation, which drives the effective tempera- 
ture to a limiting value of Xoo = T^/Tz- While Langer 
assumed Too = T, 
claim that 

The second term describes the cooling, which depends 
on the thermal frequency p. k is a free parameter that 
will give us some control of the interesting stress scales, 
as will be explained shortly. The factor exp(— /3/x) is 
the population of sites at which the effective tempera- 
ture configurations can relax, /? being the typical energy 
of such a fluctuation. These need not be STZs, but sim- 
ply any high energy conflguration, bistable or not. The 
distinction is small, and we expect /3 ^ 1, but it is nec- 
essary because steady-state localization only can occur 
when 13 ^ \. 

The A dynamics can be ignored because the STZ den- 
sity equilibrates with the effective temperature much 
faster than x equilibrates to its flnal value. The typi- 
cal timescale for A is toF"^ while for x the timescale is 
To (eoe^^/-^F) . If we assume that when T^ff = Too 
the entire material is fully covered with STZs, then 
£0 ~ ^00 ~ 6^/-^°° , and the ratio of these two timescales 
is exp^— ^ (Sa^^)^ For the metallic glasses, Tz ^ 

25,000 K, Too ~ 900 K, and the temperatures of interest 
in are all below 700 K, giving a ratio on the order of 10~^. 
The STZ dynamics are at a much faster timescale than 
the effective temperature and we are justiflcd in enforcing 
A = exp(-l/x). 
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The solutions to x = vs. stress for (3 = 0.6 for three 
different temperatures are shown in the top graph of 
Fig. 1. For certain ranges of stresses, there are three 
different effective temperatures are steady-states, which 
means a uniform system can have three possible strain 
rates. I expect the lowest and highest strain rates to be 
stable, but the middle rate unstable. A system feeling 
a stress in this unstable range, and at a strain rate be- 
tween these highest and the lowest strain rates, can be 
expected to break up into regions of fast and slow strain- 
ing. A quantitative examination is the aim of the next 
section. 



first order, the perturbation is governed by 



where Fq = r(so,'77o). Since ^ 



dt 



0, Fo 



Kpe 



-(/3-l)/X0 ( Xo-Tl 



^ I ■ Using that to remove any ex- 
plicit dependence on the stress in the stability exponent 
Lu{k) yields 



III. SHEAR LOCALIZATION: LINEAR 
STABILITY 



A. Theoretical Results 

As evidenced by Figs. 10 and 11 in [1], the STZ equa- 
tions Eqs. (10) and (17) can. under certain conditions, 
have multiple steady-state solutions corresponding to dif- 
fering effective temperatures and strain rates at same 
stress. I now examine the model to find these neces- 
sary conditions on the parameters and state variables of 
a strip of material under simple controUcd-shear-stress 
conditions such that it is unstable. Denoting the flow di- 
rection as X, by symmetry the state variables s, to, and 
X will be functions of y only. We shall assume inertial 
effects are negligible, which will allow us to assume the 
stress throughout the strip to be uniform and equal to 
the applied stress. Likewise, since the to steady state de- 
pends solely on the stress, it too will be uniform. There- 
fore any shear banding will be caused by inhomogeneities 
in the effective temperature. Consequently we shall study 
the stability of perturbations of a uniform steady state 
(so, "^0, and xo) of only the effective temperature. With 
diffusion effects included into Eq. (17) the effective tem- 
perature obeys 



X = \ e ^/'^F(s, to) (xoo - x) 



To Co 



-Pe-^/M^^-x)+^x0i. (18) 



Note that we only are considering diffusion perpendicular 
to the flow. For particles in a simple shear flow, the 
diffusion along the flow is enhanced by Taylor dispersion 
effects which are beyond the scope of this paper. 

Consider a perturbation the effective temperature of 
form X = Xo + Xi^^'^^e"''"'^*, of a system under controlled 
uniform stress sq, and with the uniform, steady-state 
state variables too and xo- To zero-th order, all of the 
time derivatives [Eqs. (10) and (18)] are equal to 0. To 



coik) = ^ 
To Co 



Kpe 



-Hlxo 



1-/3 



Xo 



T 



T 
Tz 



Xo 



Xo 



(20) 



The most unstable mode is a uniform effective temper- 
ature change (fc = 0). Generally, the effective tempera- 
ture will be between T^o and T, which, in the case where 
Too > r, means that /3 > 1 is required for instability. 
More precisely, there will only be a region of instability 
for material parameters satisfying 



(3> 



Xoo T 



(21) 



If this condition is satisfied, there will be a region of 
instability for effective temperatures between x± given 

by 



1 



l + TfTJ^X 



Ax 



(22) 



where X=\ {xoo + ^) and Ax = | (xoo - 

Interestingly, the bistability conditions arc indepen- 
dent of K. Varying the relative strength of the cooling n 
does not change the unstable effective temperature range, 
but does vary the corresponding stresses; a stronger cool- 
ing is simply overcome at a higher stress. 

While I have shown there are instabilities, I have not 
proven they are where the effective temperature as a 
function of stress is multivalued. For the a > 1, Fo is 
a monotonically increasing function of the applied stress 
So, i-C. ^^f- > 0. Therefore, the section of the curve 
where L |Ea < q requires < 0. Setting the 
LHS of Eq. (17) to and differentiating with respect to 
Xo yields 



9Fo 1-13 T 
dxo Xo V Tz 



Xo 



T 

Tz. 



Xoo - Xo 



(xw(0). (23) 
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FIG. 1: Effective temperature x a-nd instability exponent 
uj{k = 0) (in units of vs. stress s for tlie steady state 

solution of the STZ equations at varying temperatures for 
/3 = 0.6, K = 10-^ n = 24500^^, Tr = 850K, To = 836K, 
Ti = 45K, Too = 9S7K 
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FIG. 2: Analytic fit of metallic glass viscosity to the curve 
predicted from the STZ model with effective temperature, 
Eq. (25), using the Cohen-Grest expression for ce{T), Eq. (13), 
with 770 = 1.6 X 10"" Pa sec, = 24,500 K, To = 836 K, 
Tr = 850 K, and Ti = 45 K. The circles and squares are ex- 
perimental results for Vitreloy 1 from [15] and [2] respectively. 



Therefore, as one would expect, the section of the steady- 
state X vs. s curve where |^ < is unstable. 

Figure 1 shows the effective temperature and instabil- 
ity for three different temperatures. A comparison of the 
graphs of x and uj{k = 0) shows that, as predicted, the 
instability occurs at those stresses where the x is multi- 
valued. 

The system is only unstable for a finite range of strain 
rates; a material deforming sufficiently fast will be ho- 
mogeneous. Experimental evidence for such a transition 
back to homogeneous flow has been seen in some metallic 
glasses[14]. When the material is driven at an unstable 
strain rate, it can break up into two or more regions of 
different strain rates. Because the dispersion relation, 
Eq. (20), does not have a peak at nonzero fc, there is 
no length scale describing the width or spacing of the 
shear bands. The only length scale is the width a of the 
boundary between the two regions 




13/xo 



T 



T 

Tz 



XO 



(24) 



B. Comparison to Experiment 

Before doing quantitative fits to the experiments, we 
must be explicit in mapping the STZ equations as writ- 
ten to the results of a 3-D uniaxial compression test, 
for which we will follow the procedure laid out in [1]. 
The stress tensor in such tests is, denoting the axis of 
compression as x, a^x — cr' and all other components 



zero. By using the second invariant of the deviatoric 
stress tensor to define the STZ equation stress s such 
that s ~ ^j-, the STZ yield stress becomes qualita- 
tively similar to the von Mises yield criterion. Likewise, 
the scalar m will correspond to the second invariant of 
an STZ bias tensor. With these conventions, Eq. (3) de- 



scribes ixx , up to a factor of y | which can be absorbed 

into Eo, while other equations remain unchanged. 

As shown in [1], at low stresses, the material will ap- 
proximate a Newtonian liquid with a measured viscosity 
r]N{T) given by 



Vn{T) = -qo-j^cxrpl — 



a{T) 



a{T) 



(25) 



where ryo 



eo/to ' 



The metallic glass viscosities are 



separated into two groups. First the high temperature 
limit of Eq. (25), miT) = f^Texpf ^^t^ 1 is used 



Tr " \ T 

to fit the higher temperature values[lo]; then those con- 
straints are used when the lower temperature values [2] 
are fit with the full expression. The best-fit parameters 
are 7?o = 1.6 x 10"" Pa sec, Tq = 836 K, = 850 K, 
and Ti = 45 K. Tz is 24,500 K which is on the order 
of the value of iil^ for an STZ of a few atoms in diame- 
ter. The fit was done more rigorously than and improves 
upon the prcvious[l]: only two of the low temperature 
viscosities are not well-explained by Eq. (25), the first 
and the fourth. Unlike [l], I will ignore any possible er- 
rors in the temperatures reported in the experiment, and 
use the reported values. 

Before fitting the uniform steady-state behavior at 
higher strain rates, it will be useful to reduce the pa- 
rameter space by estimating some of the parameters be- 
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forehand. If we assume Too = a-s assumed previously, 
then to match the experimental observation of localiza- 
tion at T = 663 K, /3 must be less than 0.47. I feci that 
is too low and instead shall use T^c = TgoUdus which is 
937 K for the metallic glass[lG], in which case /? < 0.63. 
I shall stay at the high end of that limit and use (3 ~ 0.6. 

Next, for /2, which is assumed to be the yield stress, 
I propose the following argument: Imagine applying a 
constant strain rate to a material initially at rest. The 
effective temperature dynamics are slow enough that the 
number of zones is essential unchanged during the ini- 
tial stress rise. The peak stress is determined by the 
stress necessary to produce the applied strain rate purely 
plastically. If the strain rate is small, this can be accom- 
plished purely from plastic creep; the stress will be below 
the yield stress and highly dependent on temperature. If 
the strain rate is large, the system will have to enter the 
regime of STZ flow; the stress will be above the yield 
stress and less dependent on temperature. This exact ef- 
fect can be seen in the metallic glass data in Fig. 4 in [2], 
with a distinct change in the temperature dependence 
of those experiments where the peak stress was above 
1500 MPa. In the 2-D simulations, this corresponds to 
ii = 500^3 MPa. 

That leaves k. Instead of trying to fit uniform steady- 
state stress vs strain rate for all eight temperatures di- 
rectly, I used the minimum strain rate at localization 
as a function of temperature first. This fit is shown in 
Fig. 3. Because of the large error bars, and the vagueness 
in deciding whether a sample is localized or not, the fit 
was done by eye to yield a value of 2.5 x 10~^. While the 
model agrees well at high temperatures, at lower temper- 
atures it underestimates the strain rate at localization by 
a factor of 100. Remember, the predictions are for the 
steady state results. Later, when the dynamics are in- 
cluded, the difference will be apparent. 

Figure 4 compares the steady-state stresses vs. strain 
rates predicted by the model to those reported in the 
experiment [2]. Compared to the fit in [1] the high tem- 
perature, low stress data part of the fit is improved, but 
the low temperature, high stress parts are worse. 



100 



0.01 



0.0001 



le-oe 





1 1 1 1 1 1 1 1 1 1 1 1 1 1 


r Stable 




r ^^^""^ 
- ^--^ 


- 


!^^^ unstable 






\ 


- 






Stable 1 


' ,X , . 1 . . . . 1 . . . . 1 . . . 


. 1 .... 1 .... 1 . ; 



560 



580 



600 620 640 

Temperature [K] 



660 



680 



FIG. 3: Stability map of the uniform steady-state solutions of 
the STZ model with effective temperature for Too ~ Tsoiidus = 
937 K, /3 = 0.6, K = 2.5 X 10"^ p. = 500^ MPa. The other 
parameters are the same as in Fig. 2. Strain rates inside 
the curve are unstable. The circles are experimental data for 
Vitreloy 1 taken from Fig. 11 in [2]. 
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FIG. 4: Steady-state stress vs. strain rate for uniform solu- 
tions of the STZ model with effective temperature using the 
same parameters as Fig. 3. The symbols are the experimental 
results for Vitreloy 1 from [2] 



IV. SHEAR LOCALIZATION: NONLINEAR 
STABILITY 

A. Strain Rate Dependence 

Now that we've found the parameter ranges with the 
best prospects for banding, I will numerically solve the 
system of equations in the linearly unstable regime and 
study the full results. I will model a material under a 
controlled constant total strain rate, first starting with 
a system at an unstable steady state solution of the dy- 
namics. Then I will examine systems starting at rest, and 
the effects the transient response has on the stability. Fi- 
nally I will vary the temperature and study the effects of 
the initial STZ density. I shall assume the elastic wave 



speed is great enough that the stress can be assumed to 
be uniform at all times. The uniform stress forces the 
STZ bias ratio m to also be uniform. The time evolution 
of the stress is determined by 

s = 2fi ^£*°' -J^j^ eP^{x)dx^ . (26) 

This equation with Eqs. (10) and (18) completely de- 
scribes the dynamics of the system. 

The boundary size of any shear bands can be esti- 
mated from Eq. (24) and an estimate of the diffusion 
constant. According to Ono et al. [17], despite being 
non-equilibrium, the diffusion constant in the direction 
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perpendicular to the flow still obeys the Stokes-Einstein 
relation D oc The zone size I is on the order of a 

few atoms; I ~ 10^'' m. Despite no longer being in a 
viscous regime, the viscosity can be approximated by the 
average viscosity, j . Since the effective temperature and 
the strain rate vary greatly between the cold side, where 
X — X-j a-nd the hot side, where x = X+j the diffusion 
constant varies from 10~^^ to 10~^^ m^/sec. I shall use 
a value at the high end of this range, D = lO^^*' m^/sec, 
which is in rough agreement with the experimentally 
measured atomic diffusion constants. [IS, 19] I shall also 
ignore any temperature dependence. Any errors in this 
approximation should not seriously affect the results; as 
this is the only length scale in the problem, all the other 
lengths can be rescaled to compensate. 

Using this value in Eq. (24) gives an estimate of the 
boundary width on the order of a few angstroms. To 
avoid the divergence from the term in brackets at x±j the 
estimate Xo = f (x+ + X-) was used. It should be noted 
that although, over the temperature of the experiment, 
the strain rate and diffusion constant vary by three orders 
of magnitude, the boundary size stays relatively constant 
as the temperature varies between 577 and 663 K. This 
can be seen from rewriting Eq. (24) (ignoring the slowly 
varying, non-exponential term in brackets) as 
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The number of zones, exp(— 1/x), is the most varying fac- 
tor as the temperature decreases, but since it affects both 
the plastic work and the Stokes-Einstein approximation 
of the diffusion constant in identical ways the boundary 
size does not vary by more than a factor of 10. 

The numerical setup is as follows: The systems start 
in a uniform state given by Sq, mo, XOj a-^d have peri- 
odic boundary conditions. All systems are 0.1 microns 
in size and divided into 1000 points. Inhomogeneities 
are introduced by adding a random noise, Ax(a^j), to 
the effective temperature at each point. These random 
variables are independent and identically distributed uni- 
formly over the range {— 0.5%xo, 0.5%xo}- The system is 
then evolved according to the equations of motion up to 
strains of e/ = 0.4. From this, a final, steady-state, sta- 
tionary profile is also computed by minimizing the time- 
derivatives of the state variables as expressed in Eqs. (10), 
(18), and (26), starting with the results of the simulation 
of the dynamics at strain Sf. The same material parame- 
ters as in Fig. 3 are used with /j, having the experimental 
value of 35 GPa, and £q = exp(l/xoo)- The remaining 
parameter cq should be of order 1 and, by comparison of 
the response of the material as a shear band forms, and 
the experimental results, a value of cq = 0.6 is used. [20] 

Figure 5 shows the stress relaxation as a function of 
strain for a system with the previously listed parameters, 
at a bath temperature of 643 K, starting at the unstable 
uniform solution of Xo = 5 (X-i- + X-) = 0.313, Sq = 0.4, 
Too = 0.98. This corresponds to an applied strain rate of 
1 sec~^. 
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FIG. 5: The stress relaxation of an STZ material as a shear 
band forms. The system is started at a uniform unstable 
point, with a small amount of noise added to the effective tem- 
perature profile, then evolves at a fixed strain rate of 1 sec~^, 
according to Eqs. (10), (18), and (26). The "X" represents 
the point at which locaUzation is defined to occur. 



Localization is determined by the Gini coefficient of the 
strain rate distribution. The Gini coefficient can be de- 
fined as the average difference between all possible pairs 
of values (including pairs where both values are chosen 
to be the same point) divided by the average and nor- 
malized to be 1 in the case of an infinite system where 
only one value is nonzero. Denoting the Gini coefficient 
by (p, the value is defined by[21] 



(^.)| 



(28) 



where N is the number of sites, and e^' is the average 
of the plastic strain rate distribution. ranges from 
for a completely uniform distribution to 1 for a system 
where all the strain is occurring at one infinitesimally 
small point. Any systems with > 0.5 will be considered 
localized. This value is arbitrary, but easily measurable, 
and in most of the results (j) is sharply increasing when 
this criterion is surpassed. 

In Fig. 5, the point where the material localizes is 
marked with an "X". Any results beyond that point 
should be taken with a grain of salt, both because the 
experiments cease at that point so there is nothing with 
which to compare, but also because the sharp sudden 
onset of a shear band could lead to fracture. 

The effective temperature profile shows five shear 
bands of high strain rate and effective temperature. 
There was no a priori reason to expect five bands; with 
no relevant length scale, the number of zones is deter- 
mined by the initial noise. Similar runs were done for 
other noise realizations, and the number of bands varied 
from 2 to 5. In comparison to experiments though, the 
shear bands computed are much thinner than those re- 
ported. In a similar metallic glass, band widths of 50 nm 
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FIG. 6: Steady-state profile of the effective temperature x and 
strain rate vs. position x of the controlled-strain-rate driven 
system of Fig. 5. The lengths are in units of the system size, 
fO"'^ m. 



were reported [6] and are possibly higher in this particu- 
lar metallic glass. [22] The dashed line is the strain rate, 
which because of its strong dependence on x varies even 
more sharply than the effective temperature. 

These results are interesting and prove an instability 
is possible, but they do not directly correspond to ex- 
periments, the experiments are not started at a uniform 
shear solution. Instead they are started with no stress 
(sq = 0) and no bias in the STZ population (mg = 0). 
The stability of the initial, transient response may differ 
both qualitatively and quantitatively from the stability of 
the uniform steady-state solution. In an examination of a 
Spaepen's free volume model[7], Huang, et al.[2.3] showed 
that although the uniform steady state was stationary, an 
instability can arise during the initial response during a 
constant-strain-rate experiment. 

To study the effects of the initial loading, four addi- 
tional simulations are run. Like the previous, a system 
is created with a width of 0.1 microns discretized into 
1000 sites. Unlike before the initial stress and bias are 
zero (sq = 0, mo = 0). The initial noise of the effective 
temperature is again 0.5%xoj but xo is equilibrated to 
the bath temperature which again is chosen to be 643 K. 
This system is then used as the starting configuration for 
four controlled strain rate tests, one at each of the strain 
rates in Fig. 2 from [2]. 

The stress as a function of strain for these simulations 
is plotted in Fig. 7. Like the metallic glass experiment, 
the lowest three systems stay homogeneous throughout 
the run, while the highest strain rate localizes suddenly. 
The graphs for the two highest of non-localizing strain 
rates, 3.2 x 10"'^ and 5.0 x 10"'^ sec~^, display a stress 
overshoot, also like the experimental data. 

The highest strain rate results also shows a stress over- 
shoot quickly followed by the onset of localization, both 
in the numerical results and in the actual experiment. 
Final steady-state effective temperature profiles arc com- 
puted by the same procedure as previously explained and 
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FIG. 7: Stress vs. strain for a system initially at rest for 
four different strain rates at T = 643 K. In addition to 



the parameters of figure 3, co = 0.6, eo 
10"^°mVs (^^) = 1.12 * 10"^°, and the number of lat- 
tice points is 1000. 
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FIG. 8: Steady-state effective-temperature profiles of the sys- 
tems in Fig. 7. 



plotted in Fig. 8. Note how the system at the higher 
strain rate can find a lower average effective temperature 
by localizing the shear. The strain rate is so much larger 
in the shear band, the remainder of the material can relax 
to a lower effective temperature. 

Unlike the system that started in the uniform steady 
state, the time it takes for the highest strain rate system 
in Fig. 7 to shear localize is much shorter, less than 1 sec- 
ond compared to the 5 seconds for a system at the most 
unstable strain rate. The speed is even more dramatic 
when considering that the dominating time scale is the 
strain rate, which is 10 times faster in the earlier case. 

A similar localization-during-transient response was 
seen by Huang, et al.[2.3] during their numerical analysis 
of Spaepen's free- volume model of plasticity [7] . Although 
that model has only stable uniform steady states, they 
showed that the transient response of such a system could 
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FIG. 9: Stability exponent lo{0) given by Eq. (20), for the 
uniform solution of the STZ equations as a function of strain 
for a material initially at rest and under an applied strain rate 
of e = 10~^s~^, the highest strain rate in Fig. 7. 

lead to a localization. (Their localization was only tem- 
porary, but that need not the case with the STZ model 
with effective temperature.) Following their logic, a time- 
dependent stability exponent can be calculated by exam- 
ining the stability of the uniform solution as the latter 
evolves in time [24]. Instead of perturbing around a time- 
independent solution, the perturbation is done with re- 
spect to a time-dependent uniform solution SQ{t), moit), 
and Xo{t)- The result is exactly the same as Eq. (19), 
but no longer can the stress-dependence be explicitly re- 
moved as in Eq. (20). Figure 9 shows the largest stability 
exponent uj{0) taken from Eq. (19) as a function of strain 
rate for a uniform system with the same parameters as in 
Fig. 7 at a strain rate of 10~^ s~^. Because of the initial 
stress rise, T becomes very large near the stress peak and 
the solution becomes very unstable to localization. In 
comparison to the final stability exponent of the uniform 
solution, the peak value is over two orders of magnitude 
larger. 



B. Temperature Dependence 

Now that we've examined the effects of varying the 
strain rate, we will vary the bath temperature as well. 
To do so, first more consideration must be made as to 
the initial values of the effective temperature. So far, 
we have examined the model for cases where the initial 
effective temperature is equal to the bath temperature. 
If this were always the cases, there would be no reason 
that the effective temperature had to be "effective" and 
couldn't actually be a traditional thermal temperature. 
In which case, all the results seen so far would simply be 
examples of adiabatic shear bands. The only difference 
would be the diffusion constant, but that only sets the 
length scale. But as the initial temperature lowers, that 
equality will no longer hold. Instead we must estimate 
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FIG. 10: Initial value of the effective temperature xo as a 
function of bath temperature for a material obeying the STZ 
equations with effective temperature with the same parame- 
ters as Fig. 5 after being cooled from TiniUai = Teff = Ti„f 
to 295 K at 1 K/s, then, some time later, heated to the final 
temperature T at 20 K/min, and held at that temperature 
for 10 minutes. 



the initial effective temperature from the thermal history 
of the material. 

Consider a material that is first cooled from the melt- 
ing temperature to room temperature, then heated to 
and annealed at a temperature at which it will be tested. 
Initially, Tf,ff = T, but as T drops, the cooling dynamics 
become slower and slower, and eventually the effective 
temperature of the system falls out of thermal equilib- 
rium at some value. Assuming a final bath temperature 
of 295 K and a cooling rate of 1 K/s (the typical cooling 
rates in the production of Vitreloy 1 are between 10 and 
0.9 K/s [16, 25]), and using the same parameters as in 
the previous section, the effective temperature plateaus 
at 640 K after roughly 20 minutes of cooling and stays 
frozen at that value for much, much longer than labora- 
tory timescales (over lO""^" years). 

Immediately before testing, the samples are heated to 
the testing temperature and annealed for a few minutes. 
During the annealing process, if the final temperature is 
high enough, the effective temperature will have time to 
equilibrate to the bath temperature. Figure 10 shows the 
final effective temperature as a function of the bath tem- 
perature for a system following a thermal history similar 
to that performed in [2], cooled from the melt to room 
temperature, held there for more than a few hours, then 
heated at a rate of 20 K/min to a final temperature at 
which they are anneled for 10 minutes. At high tem- 
peratures, the system completely equilibrates, and the 
thermal history is forgotten; at low temperatures, the 
thermal effects are still not strong enough to noticeably 
change the temperature; and in between the material is 
equilibrating over laboratory timescales. 

Figure 11 plots the stress as a function of strain for a 
strain rate of 0.1 at five different temperatures with 
the same parameters as used in Fig. 7, except for the 
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FIG. 11: Stress vs. strain for a system initially at rest for five 
different temperatures at a strain rate of e = 0.1 s~^. All 
other parameters are the same as in Fig. 7. 



initial effective temperature, which is taken from Fig. 10. 
Note that the initial heterogeneities are calculated as a 
percentage of the initial effective temperature; if a point 
is 0.25% above the average at T = 295 K, it is also 0.25% 
above the average for T = 683 K. 

Qualitatively, the model agrees well with experiment. 
Localization is only seen in at the lowest three tem- 
peratures, occurring when the stress sharply drops. 
Quantitatively, the peaks are higher than seen in the 
experiments [2], 3150 MPa from the model versus 1850 
MPa in the experiments, but this is the regime of low 
temperatures and high strain rates where the rate fac- 
tors are most likely to be incorrect. This though is much 
better agreement than if the effective temperature was 
taken to be the actual bath temperature. In that case, 
the initial number of zones for the lowest two temper- 
atures, the stresses become outrageously bigger than in 
the experiments. Roughly 5,000 MPa for T = 523 K and 
25,000 MPa for T = 295 K. Trying to ehminate the effec- 
tive temperature by simply using the bath temperature 
will not work. 

Figure 12 shows the peak stress vs. strain rate for 
a number of different temperatures. Solid symbols are 
those that remain localized up to the final strain; open 
symbols exhibit strain localization. The logic used earlier 
predicts that those runs where the peak stress goes above 
1 should localize and the numerical results support this. 
All the runs with peak stress below 0.974 stay uniform; 
all the runs with peak stress above 0.99 localize. (There 
were no results in between.) Any samples whose peak 
stress rises above this value localize, while any materials 
where the stress stays below it remain uniform. When 
the strain rate is low and the temperature is high, the 
response is dominated by linear viscosity, as expected. 
For the higher stresses, the material is forced at a rate 
higher than the creep alone can deliver. Instead, the flow 
is dominated by the STZ flow, and the temperature be- 
comes irrelevant. More important is the initial number 
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FIG. 12: Peak stress vs. applied strain rate at various tem- 
peratures for STZ models initially at rest. The dashed line is 
the STZ yield stress. 



" 3000 

■i 

g_2000 







1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 












G— O j 


= 0.0001 s"' 






□ — □ 


0.00032 






0—0 


0.0032 








0.01 








0.1 








200 










*s — < — *>sl^ 




O 








y 






1 




. . 1 , 


, , , 1 .... 1 . VA 



300 400 500 600 700 

temperature [K] 

FIG. 13: Peak stress vs. temperature at various applied strain 
rates for STZ models initially at rest. The dashed line is the 
STZ yield stress. 



of zones. 

Figure 13 shows the peak stress as a function of tem- 
perature for varying strain rates. Again, the change in 
behavior as the peak stresses rises above the yield stress 
is seen. Note that although the data for e = 200 in 
Fig. 13 appears to become viscous at high temperatures, 
the change in slope is a result of the change in xq. Re- 
member, above 620 K, the initial effective temperature is 
equal to the bath temperature, so when the temperature 
increases, not only does p increase, but xq does as well. 
More zones are available to contribute to the plastic flow, 
meaning the same strain rate can be attained at a lower 
stress. 

Figure 14 is a scatter plot of the peak stress at many 
different pairs of temperature and applied strain rate. 
The circles are the experimental strain rates below which 
the system is uniform and above which the system local- 
izes taken from [2]. Inside the solid line is the region in 
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FIG. 14: Scatter plot of temperature and applied strain rate 
pairs for the STZ model starting at rest (so = 0, mo = and 
Xo according to Fig. 10). The open squares correspond to runs 
that exhibit localization while the solid squares remained uni- 
form. The circles are experimental data. The solid line is the 
boundary of instability under infinitely small perturbation. 
The dashed line is emaxcreep{T). Below it, the flow is pri- 
marily creep, and above the flow is primarily from the STZ 
dynamics. 

which the steady state of the model is non-uniform. The 
solid squares correspond to systems in which the strain 
rate remains uniform up to a strain of 0.45; the open 
squares are those systems where the deformation local- 
izes. The data shows that the existence of a shear band is 
highly dependent on the initial conditions. For example, 
while all the points tested below 500 K should be uniform 
in the steady state , they all localize as the system re- 
sponds to the applied strain rate. Instead of a a uniform 
increase in effective temperature towards the final uni- 
form solution, the material creates one large shear band 
first where the effective temperature is greatest. That 
shear band allows the stress to drop, slowing the heat- 
ing in the remainder of the system. Instead, the system 
is heated by diffusion of effective temperature from the 
shear band and the band widens to absorb the entire sam- 
ple. The numerical equations appeared too stiff to verify 
the eventual return to homogeneous flow at that temper- 
ature, but it was verified for T — 690 K and e — 200 sec~^ 
where (j) drops below 0.5 at a strain of roughly 15. Like- 
wise, the point T = 593 K and i = 3.2 x 10~^ sec~^ which 
was initially uniform by my criteria, slowly but steadily 
localizes, reaching = 0.5 at a strain of 0.9. 

If the logic proposed earlier holds true, those strain 
rates that are low enough that they can be produced 
solely by the thermal creep from just the initial density 
of zones should not show localization in experiment. In 
other words, for i such that 

£ < ema.creep{T) = ££ ^Z^^" ^^^C (1) {T{1} - m,,(l)) , 
To 

(29) 

the material should appear to remain uniform, while ma- 
terials forced at strain rates i > £maxcreep{T) should lo- 



calize. £maxcreep{T) is the dottcd liuc and it agrees well 
with the numerical results, only differing for those points 
where the max stress is between 0.974 and 0.99 s„. 



V. CONCLUSION 

As shown in the previous sections, the addition of the 
effective temperature dynamics into the STZ theory en- 
ables it to qualitatively describe the localization phe- 
nomenology seen in experimental data. 

The shear banding caused by this characteristic of the 
model is different from that proposed in [26] . In that pa- 
per, the bands were the result of the interactions between 
STZ flips mediated by the elasticity of the material. The 
resulting shear bands were characterized by differences in 
the STZ bias, A = mA. In contrast, these shear bands 
are caused solely by inhomogeneities in the density, not 
bias of STZs. The previously proposed mechanism leads 
to bands of a definite wavelength, while the mechanism 
discussed here docs not. The wavelength undetermined, 
but the number of bands itself can vary depending on the 
initial conditions. 

In addition to an instability of the steady state, the 
transient response to initial loading can develop shear 
bands. When a system where the effective temperature 
is noisy is loaded at a constant strain rate, the "hot- 
ter" regions, because they have more zones, experience a 
disproportionate share of the plastic strain. More defor- 
mation means more dissipation from plastic work, which 
means the effective temperature increases faster, which 
means more zones, which allows for greater plastic defor- 
mation. As a result, the deformation can possibly localize 
in some regions, while the remainder of the material is 
relatively undeformed. 

This instability in the response of the material to the 
initial loading was shown to well describe the strain rates 
and temperatures for which localization appears in the 
experiments. Since this dynamical instability map does 
not necessarily require instabilities of uniform steady- 
state results, such an instability might exist in the origi- 
nal /3 = 1 version of the model Langer studied previously. 
That version does fit the steady-state stress vs. strain rate 
data better and is a simpler form in which to compute 
results analytically. Such a study is an obvious avenue 
for future research. 
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